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The coupled-channels technique is applied to analyze recent fusion data for 48 Ca+ 48 Ca. The 
calculations include the excitations of the low-lying 2 + , 3~ and 5~ states in projectile and target, 
and the influence of mutual excitations as well as the two-phonon quadrupole excitations is also 
investigated. The ion-ion potential is obtained by double-folding the nuclear densities of the reacting 
nuclei with the M3Y+repulsion effective interaction but a standard Woods-Saxon potential is also 
applied. The data exhibit a strong hindrance at low energy compared to calculations that are 
based on a standard Woods-Saxon potential but they can be reproduced quite well by applying 
the M3Y+repulsion potential with an adjusted radius of the nuclear density. The influence of the 
polarization of high-lying states on the extracted radius is discussed. 

PACS numbers: 24.10. Eq, 25.60.Pj,25.70.-z 



I. INTRODUCTION 

Heavy-ion fusion reactions are sensitive probes of the 
nuclear surface of the reacting nuclei. Roughly speaking, 
the height of the Coulomb barrier is determined by the 
radii and diffuseness of the densities, and the enhance- 
ment of fusion at subbarrier energies is governed by cou- 
plings to the excitation of low-lying surface modes [![. 
This picture may not always succeed in reproducing the 
measured cross sections. In practical coupled-channels 
calculations it is often necessary to make adjustments, 
either in the structure input or in the ion-ion potential. 
The adjustments may reflect the influence of high-lying 
states or other reaction channels that are not treated 
explicitly in the calculations. Another complication is 
the hindrance of fusion which occurs at low energies and 
very small cross sections [2\. The hindrance can be ex- 
plained, for example, by using a shallow potential in en- 
trance channel [3j or by modeling the fusion dynamics 
for touching and overlapping nuclei [J] . 

In this work the coupled-channels technique is applied 
to analyze the fusion data for 48 Ca+ 48 Ca that were re- 
cently measured down to very small cross sections below 
1 /do [Sj. An analysis of the data provides the opportu- 
nity to investigate whether the fusion hindrance, which is 
a well established phenomenon in extreme subbarrier fu- 
sion reactions of medium-heavy systems with large neg- 
ative Q-values @, also occurs in the fusion of calcium 
isotopes with near zero Q-values. An indication of a hin- 
drance in the fusion of 48 Ca+ 48 Ca has already been ob- 
served [1[ because the low-energy data could only be re- 
produced by coupled-channels calculations that employ a 
very large diffuseness of the ion-ion potential. However, 
the hindrance observed there was not strong enough to 
show an S factor maximum in the energy region of the 
measurement. 

The ion-ion potential and the nuclear couplings that 
will be used are derived as in previous work [3[ from the 
double-folding of the densities of the reacting nuclei and 
the M3Y+repulsion effective interaction. Since the low- 



lying structure of 48 Ca is fairly well established and the 
influence of transfer reactions is always suppressed for 
symmetric systems, it is expected that the simple picture 
of fusion described above would apply to the fusion of 
48 Ca+ 48 Ca. On the other hand, it is well known that the 
excitation and/or polarization of high- lying states that 
are included explicitly in the coupled-channels calcula- 
tions can lead to a negative energy shift of the calculated 
fusion cross sections [f|. This is effectively equivalent to 
increasing the radii of the reacting nuclei. The radius 
that is extracted by optimizing the fit to the fusion data 
can therefore be too large because it can be contaminated 
by the polarization of states which are not included ex- 
plicitly in the calculations. It is of interest to see how the 
extracted radius depends on the model space of excited 
states that are included in the calculations, and how well 
it compares to the expected matter radius of 48 Ca. 

The study of the fusion of different calcium isotopes 
started more than twenty years ago with the measure- 
ments by Aljuwair et al. [7| but the cross sections were 
only measured down to about 1 mb. The most challeng- 
ing theoretical issue at that time was to explain the fusion 
of 40 Ca+ 48 Ca which appeared to be strongly enhanced 
by couplings to transfer^reactions, in particular to those 
with positive Q values [1, Q . A strong motivation for re- 
viving the study of the fusion of calcium isotopes is that, 
in addition to 48 Ca+ 48 Ca [|, the fusion of 40 Ca+ 48 Ca 



has also recently been measured down to the 1 /jb 1C | , 



and a new measurement for 40 Ca+ Ca is underway 
In order to be able to focus on and isolate the effect of 
transfer on the fusion of the asymmetric 40 Ca+ 48 C sys- 
tem, it is necessary first to develop a good description of 
the fusion of the two symmetric systems, and the present 
work is a step in that direction. 

The nuclear structure properties of 48 Ca are discussed 
in the next section and Sect. Ill describes the construc- 
tion of the ion-ion potential. The coupled-channels tech- 
nique is summarized in Sect. IV together with an analy- 
sis of the data that is based on Woods-Saxon potentials. 
The analysis based on the M3Y+repulsion potential is 
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presented in Sect, 
in Sect. VI. 



V. Finally, the conclusions are given 



II. NUCLEAR STRUCTURE INPUT 



The nuclear structure input to the coupled-channels 
calculations is shown in Table HI The elastic channel and 
the one-phonon excitations of the low- lying 2 + , 3 _ and 
5~ states in projectile and target results in a total of 
7 channels.. The coupling stren gths for the excitation of 
these states are taken from Ref. [12] where they were cal- 
ibrated by analyzing the elastic and inelastic scattering 
of 16 O on calcium isotopes [l3j. It should be noted that 
the Coulomb and nuclear coupling strengths, expressed 
in Table U by the values of f3R/\/4Tr, are different. The 
Coulomb couplings are in most cases consistent with the 
currently adopted electromagnetic transition probabili- 
ties or B-values [l4{ that are quoted in the third column 
of Table U 

Also shown in Table U are the + , 2 + and 4 + mem- 
bers of two-phonon quadrupole excitation. The adopted 
B(£'2)-values [14J shown in the third column of the sec- 
ond part of Table fl] can be combined into an effective 
two-phonon excitation. For example, the effective B(E2) 
value for the two- to one-phonon transition is given by 
the sum, 



B(E2, 2ph -> lph) = (2020|/0) 2 B(E2,I -> 2), 

7=0.2,4 

and the two-phonon excitation energy is obtained as the 
weighted average of the individual two-phonon excita- 
tions energies [151 ] . The parameters obtained for the effec- 
tive two-phonon quadrupole excitation are shown in the 
last line of Table |U Including the effective two-phonon 
quadrupole excitations in the coupled-channels calcula- 
tions, in addition to the 7 channels mentioned above, 
leads to a total of 9 channels. Unfortunately, nothing is 
known about the two-phonon excitations of the 3~ and 
5~ states so they will be ignored. 

To get a feeling of the influence of higher-order excita- 
tions on the calculated fusion cross sections one can also 
include all of the 15 mutual excitations channels that are 
generated from the six one-phonon excitations presented 
in Table |U Together with the basic 9 channels mentioned 
above, that sums up to a total of 24 channels. This will 
be referred to as the full calculation and the results will 
be compared to the fusion data and to calculations that 
include the 9 channels described above, as well as the 
no-coupling limit in which case there is only 1 channel. 



TABLE I: Nuclear structure input for 48 Ca. Values marked 
with * are from Ref. [j^| . The B(E\) values are from Ref, 
[3]. They are consistent with the Coulomb excitation pa- 
rameters (/3-R)c/v / 47r, except for the 3~ state where Ref. [12] 
uses a larger value, B(E3) = 6.8 W.u. The second part of 
the Table shows the 3 transitions that determine the effec- 
tive two-phonon quadrupole state. Its excitation energy and 
couplings to the 2f state are shown in the last line. 
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B(EX) (W.u.) 




(/3 ^_~ (fin) 

V47T 
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3.832 


1.71(9) 


0.126* 


0.190* 


3" 


4.507 


5.0(8) 


0.250* 


0.190* 


5" 


5.146 




0.049* 


0.038* 


0+^2+ 


4.283 


10.1(6) 


[0.098] 




4+^2+ 


4.503 


0.261(6) 


[0.025] 




2+^2+ 


5.311 


9(9) 


[0.111] 




Eff 2PH 


4.849 


4.7(29) 


0.15 


0.15 



III. THE ION-ION POTENTIAL 

The parameters of the Woods-Saxon (WS) potentials 
that will be used in this work, 



Uws(r) 



■V 



l + exp((r-i? )/a)' 



(2) 



are those proposed in Ref. Eqs. (III.2.40-45). The 
parameters for the system 48 Ca+ 48 Ca are a = 0.662 fm 
for the diffuseness and Vq — 64.10 MeV for the depth of 
the potential. We refer to this potential as the 'standard' 
WS potential because its diffuseness is consistent with the 
analysis of elastic scattering data (l6| . The radius i?o, on 
the other hand, will be treated as a free parameter and 
it will be adjusted to optimize the fit to the data in the 
coupled-channels calculations that are discussed in the 
next section. 

There is an interesting point concerning the isospin 
dependence of the nuclear potential proposed in Ref. [l6| . 
It enters through the average nuclear surface tension 7 
of the reacting nuclei, 



7 = 0.95 



1 - 1. 



N„ 



Z a N b -Z b 



A a 



A b 



MeV fm" 



(3) 



according to Ref. [Tg, Eq. (III. 2. 30). The correction 
factor in Eq. ©, which depends on neutron excess, is 
equal to one in reactions that involve 40 Ca. In the case 
of 48 Ca+ 48 Ca, the correction factor reduces the surface 
tension by 5%; this correction is included in the value of 
the depth parameter Vq mentioned above. 

The M3Y+repulsion potential is calculated using the 
double- folding technique described in Ref. [3]. It is 
based on the Reid parametrization of the M3Y effective 
nucleon-nulceon interaction [17J. The spherical nuclear 
densities are parametrized by 



p(r) 



p exp(i?/a) 



cosh(r/a) + cosh(i?/a) 



(4) 
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where R and a are the radius and diffuseness parame- 
ters, respectively, and po is a normalization constant. It 
is seen that the density, Eq. (j4|), approaches the Fermi 
function density po/(l + exp((r — R)/a)) for R/a >> 1. 
The parametrization Eq. (|4]) is used here because it has 
some very useful analytic properties as discussed in the 
Appendix of Ref. [llj]. For example, the Fourier trans- 
form has an analytic form which simplifies the calculation 
of the double-folding potential in Fourier space from the 
expression Q, 

U{r) = ^2 J k 2 dk p(k) p(k) v nn (k) j (fcr), (5) 

where v nn (k) is the Fourier transform of the effective 
nucleon-nucleon interaction, and jo(x) — sia(x)/x is 
a spherical Bessel function. Another advantage of the 
parametrization Eq. (j4]) is that the mean-square radius 
is given by the simple expression 

(r 2 ) = l(R*+ 7 -(anf). (6) 

The repulsive part of the M3Y+repulsion potential is 
determined by two parameters, namely, the strength v r 
of the contact effective interaction, 

v^(r)=v r 6(r), (7) 

that generates it, and the diffuseness a r of the densities 
that are applied in the double- folding calculation Q . The 
radius parameter R of the densities, on the other hand, 
is kept the same as in the calculation of the direct and 
the exchange parts of the M3Y double-folding potential. 
The diffuseness of the density that is used in calculating 
the direct and the exchange part of the M3Y potential is 
kept fixed with the value a = 0.54 fm. 

The two parameters a r and v r of the repulsive part of 
the potential are constrained so that total nuclear poten- 
tial energy, Un(t), for completely overlapping nuclei is 
consistent with the equation-of-state. That leads to the 
relation Q, 

U N (r = 0)^^K, (8) 

where A p is the mass number of the smaller nucleus and 
K is the nuclear incompressibility. For 48 Ca+ 48 Ca the 
value K = 223.7 MeV predicted by the Thomas-Fermi 
model of Myers and Swiatecki [19[ will be used. Thus 
there are essentially only two free parameters of the 
M3Y+repulsion interaction, namely, the radius R and 
the diffuseness parameter a r . They will be adjusted to 
optimize the fit to the fusion data. The strength v r of the 
repulsive interaction, on the other hand, is constrained 
for given values of R and a r by the nuclear incompress- 
ibility according to Eq. (JSJ. 

Some of the entrance channel potentials that are used 
in this work are illustrated in Fig. [TJ The height of the 
Coulomb barrier is essentially the same for all four poten- 
tials but the thickness of the barrier is very different. The 
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FIG. 1: (Color online) Entrance channel potentials for 
48 Ca+ 48 Ca obtained from the Woods-Saxon (WS, with R 
= 8.562 fm) and the pure M3Y (M3Y-1) potentials. The 
two upper curves are the shallow M3Y+repulsion potentials 
determined in Sect. V. The energy of ground state of the 
compound nucleus 96 Zr is also indicated. 

(blue) dashed curve is the entrance channel potential for 
the Woods-Saxon potential. It has the minimum pocket 
energy V m in =23.09 MeV and the height of the Coulomb 
barrier is Vcb = 51.63 MeV. The latter potential was 
determined by optimizing the fit to the fusion data with 
center-of-mass energy larger than 50 MeV. This energy 
cut was chosen because the fusion hindrance phenomenon 
sets in below 50 MeV as we shall see in the next section. 

The upper two curves in Fig. [TJare the M3Y-|-repulsion 
entrance channel potentials that are obtained in Sect. V. 
They were determined by optimizing the fit to the fusion 
data in coupled-channels calculations that include the 24 
channels described in Sect. II. There are two solutions, 
the M3Y+rep-l and M3Y+rep-2 potentials, which are 
discussed in detail in Sect. V. These potentials are shal- 
lower than the standard Woods-Saxon potential, which is 
a characteristic feature of the M3Y+repulsion potentials 
that have been extracted from fusion data 0|. 

Finally, the entrance channel potential for the pure 
M3Y(+exchange) potential is also shown. It is unrealistic 
because it is deeper than the ground state energy of the 
compound nucleus 96 Zr which is indicated by the thick 
horizontal line. 



IV. COUPLED-CHANNELS CALCULATIONS 

The coupled-channels calculations are performed in the 
rotating frame approximation, and the fusion cross sec- 
tions are determined by imposing ingoing- wave boundary 
conditions at the position of the minimum of the pocket 
in the entrance potential. This procedure is commonly 
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used and is described, for example, in Refs. [1, In 
the present work no imaginary potential will be applied. 
The fusion cross section will therefore vanish when the 
center-of-mass energies is lower than the minimum en- 
ergy of the pocket in the entrance channel potential. 

The nuclear potential enters the coupled equations 
both directly by determining the entrance channel poten- 
tial and indirectly by determining the nuclear couplings 
to first and second order in the deformation amplitudes 
through the first and second derivatives of the nuclear 
potential (see Ref. [20] for details.) There are in prin- 
ciple couplings of even higher order and to higher-lying 
states [2l[ but they will be ignored in the present study, 
partly because they are poorly known and partly because 
they are not expected to play a large role in the fusion of 
the not-so-heavy system 48 Ca+ 48 Ca. This expectation is 
based on the experience gained in Ref. [l5|. However, the 
polarization of high- lying states @ that are not included 
in the calculations could distort the analysis. 

The fusion data for 48 Ca+ 48 Ca Q are compared in Fig. 
[5] to two coupled-channels calculations that are based on 
standard Woods-Saxon potentials [l6| with the diffuse- 
ness a = 0.662 fm and depth V = 64.10 MeV. All 24 
channels described in Sect. II were included in the cal- 
culations. It is seen that the data are hindered at low 
energies and the energy dependence is much steeper than 
predicted by the calculations. The dashed curve in Fig. 
[2] is the best fit to all data points; it is achieved with 
radius i?o = 8.495 fm but the fit is very poor with an 
average % 2 per data point of x 2 /N — 7.9, including the 
statistical uncertainties and a systematic error of 7%. 

The solid curve in Fig. [5] is based on the slightly larger 
radius, Rq = 8.562 fm. It provides a better account of the 
data near and above the the Coulomb barrier as discussed 
below. The associated entrance channel potential is the 
(blue) dashed curve shown in Fig. [TJ 

The behavior of the hindrance in the fusion of 
48 Ca+ 48 Ca is illustrated in Fig. [3] in terms of the ra- 
tio of the measured and calculated fusion cross sections. 
It is seen that the ratio with respect to the best fit to the 
data (the solid diamonds) has a strong peak near 50 MeV, 
slightly below the Coulomb barrier which is at Vcb ~ 52 
MeV. The ratio drops quickly at energies below 50 MeV. 
This is attributed to the fusion hindrance phenomenon. 
In fact, the steep falloff with decreasing energy observed 
in the comparison to standard coupled-channels calcu- 
lations was the signature that was first used to identify 
the fusion hindrance [22j . Later it was shown that the 
hindrance is often so strong that the S factor for fusion 
develops a maximum at very low energies. Moreover, it 
was realized that an S factor maximum together with the 
energy Eg of the maximum is a good quantitative way 
to characterize the fusion hindrance phenomenon [23| . 

Since the fusion hindrance occurs at low energies one 
may exclude the low energy region and focus on repro- 
ducing the data at higher energies. The result of this 
approach is shown in Fig. [5] by the solid curve which 
is based on a slightly larger radius, Rq = 8.562 fm. 
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FIG. 2: (Color online) Fusion cross sections for 48 Ca+ 48 Ca. 
The error bars reflect the statistical uncertainties. The curves 
are coupled-channels calculations with 24 channels that use 
Woods-Saxon (WS) potentials with two different radii. 
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FIG. 3: (Color online) Ratios of the measured and calculated 
cross sections shown in Fig. [2] The error bars on the open 
diamonds were determined by the statistical uncertainties and 
a 7% systematic error. 



The larger radius implies larger cross sections below the 
Coulomb barrier but that gives a better description of the 
excitation function in the barrier region. The radius of 
the Woods-Saxon potential was therefore chosen so that 
the ratio of the measured and calculated cross sections 
essentially is a constant above 50 MeV. This is illustrated 
by the open diamonds Fig. [3] It is seen that the fusion 
hindrance sets in very strongly below 50 MeV, where the 
ratio falls off very steeply with decreasing energy. 
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FIG. 4: (Color online) Results of the \ 2 analysis of the 
48 Ca+ 48 Ca fusion data 0]. The x'/N, minimized with re- 
spect to the radius R, is shown as function of the diffuseness 
parameter a r . The solid curve (24ch) is for calculations that 
include 24 channels; the dashed curve (9ch) is for 9 channels. 



V. ANALYSIS BASED ON THE 
M3Y+REPULSION POTENTIAL 

The parameters of the M3Y+repulsion potential that 
provides the best fit to the data were determined using 
an improved calibration procedure. For a given nuclear 
radius parameter R of 48 Ca and a given diffuseness a r 
of the density used in calculating the repulsive part of 
the potential, the strength of the repulsive term v r was 
adjusted so that the incompressibility K = 223.7 MeV 
was achieved in Eq. (|5J|. Having determined the nuclear 
potential, coupled-channels calculations were performed 
and the average x 2 per data point, x 2 /N , was calculated 
from the statistical uncertainties and a systematic error 
of 7%. 

The above procedure was repeated with different val- 
ues of the radius R for a fixed diffuseness parameter a r 
until a minimum /N was found. The whole process was 
repeated for a new value of a r . The results of this process 
are illustrated in Fig. @]where the x 2 /N , minimized with 
respect to the radius R, is plotted as a function of the dif- 
fuseness parameter a r . The dashed curve is the result of 
calculations that include the 9 channels. The solid curve 
is the result obtained with all 24 channels described in 
Sect. II. 

The fusion data of Ref. @ are compared in Fig. [5] 
to various calculations. The solid curve is the result of 
coupled-channels calculations associated with the deep- 
est minimum in Fig. @] which has a x 2 P er data point of 
1.66. The upper (blue) dashed curve is the cross sections 
obtained in similar calculations using the Woods-Saxon 
potential with the radius Rq = 8.562 fm. All 24 chan- 
nels described in Sect. II were included in both sets of 
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FIG. 5: (Color online) Fusion cross sections for 48 Ca+ 48 Ca. 
The upper two curves are coupled-channels calculations with 
24 channels that are based on the Woods-Saxon (WS) and 
the M3Y+repulsion (M3Yl+rep-l and M3Y+rep-2) poten- 
tials discussed in the text. The lowest dashed curve is the 
no-coupling limit based on the same M3Y+repulsion poten- 
tial. 

calculations. The only difference between the two cal- 
culations is the choice of the nuclear potential, and it is 
seen that the shallow M3Y+repulsion potential labeled 
M3Y+rep-2 is a much better choice. 

The lowest dashed curve in Fig. [5] is the cross sec- 
tion obtained in the no-coupling limit (i. e., with only 1 
channel) using the same M3Y+repulsion potential that 
was used to produce the solid curve. Comparing the two 
curves, it is seen that the effect of the couplings to the 
24 channels is equivalent to shifting the no-coupling limit 
almost 1 MeV to lower energies. 

A. Details of the analysis 

The minima of the curves shown in Fig. U define the 
stable solutions of the x 2 analysis of the fusion data since 
they are minima with respect to variations in both a r and 
R. There are two local minima for each set of calculations 
and the parameters of the M3Y+repulsion interactions 
that determine them are given in Table HU It is seen the 
two solutions obtained with 9 coupled have almost the 
same x 2 /N ■ It is not clear what causes the existence of 
two solutions. The main difference between them is that 
the energy of pocket in the entrance channel potential, 
Vmin, is about 8 MeV deeper in the solution with the 
smaller radius R. 

Of the two solutions obtained with 24 channels, the 
one with the larger radius gives a much better fit to the 
data with x 2 /N — 1.66 and the associated potential will 
be referred to as the M3Y+rep-2 potential. The poten- 
tial for the solution with the smaller radius is called the 
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TABLE II: Parameters of the M3Y+repulsion potential asso- 
ciated with the x 2 minima in Fig. [?] Results are shown for 
coupled-channels calculations that include 9 and 24 channels, 
respectively. The last three columns show the minimum of 
the pocket Vmin, the height of the Coulomb barrier Vcb, and 
the \ 2 P er data point. 



No. of Ch. 


7? 


a r 


v r 




Vcb 


X 2 /N 




(fm) 


(fm) 


(MeV fm 3 ) 


(MeV) 


(MeV) 




9 


3.775 


0.4025 


480.1 


33.58 


51.67 


2.76 


9 


3.810 


0.4250 


504.2 


41.66 


51.60 


2.69 


24 


3.745 


0.4070 


481.8 


34.61 


51.77 


2.52 


24 


3.798 


0.4295 


505.6 


42.55 


51.73 


1.66 



TABLE III: The radius R of Ca extracted from the analysis 
of the fusion data; c. f. Table [Hi The rms radii are compared 
to the measured rms charge [24| . proton [24[ , and neutron [SfJ 
radii. The latter two have been combined into the rms matter 
radius shown in the last line. The quoted matter-radius ([3.75] 
fm) was derived by inserting the rms matter-radius and the 
diffuseness a — 0.54 fm into Eq. ((6J. 



Reference 


No of Ch. 


R (fm) 


{r 2 ) 1/2 (fm) 




9 


3.755 


3.547 




9 


3.810 


3.569 


M3Y+rep-l 


24 


3.745 


3.528 


M3Y+rep-2 


24 


3.798 


3.562 


Charge [24] 






3.474(1) 


protons 






3.387(1) 


neutrons [251 






3.63(5) 


matter 




[3.75] 


3.53(3) 



M3Y+rep-l potential. The two entrance channel poten- 
tials are illustrated in Fig. 1. An important question is 
whether the parameters of the stable solutions are real- 
istic, or whether some of them can be ruled out as being 
unrealistic. One parameter of particular interest is the 
radius which is examined below. 

The rms (root-mean-square) radii obtained for the sta- 
ble solutions are shown in the fourth column of Table 
IIIII They can be compared to the estimated experimen- 
tal rms matter radius of 48 Ca shown in the last line of 
the Table. The estimate was based on the rms radius 
of the proton distribution, which was obtained from the 
measured rms charge radius [24| . and the experimental 
rms radius of the neutron distribution [25j . The neutron 
radius is uncertain and several values exist in the litera- 
ture. The experimental value chosen here was obtained 
from an analysis of elastic proton scattering data at 800 
MeV [25[ and is in fairly good agreement with most of 
the theoretical predictions shown in Table 1 of Ref . (2|| . 

The estimated rms matter radius quoted in the last 
line of Table IIIII is in perfect agreement with the rms 
radius associated with the M3Y+rep-l solution. The rms 
radius for the M3Y+rep-2 solution is larger but it is still 
consistent with the experimental estimate within the lex 



TABLE IV: The \ 2 / N (column 3) obtained in calculations 
that use the M3Y+rep-2 potential and include 1, 9 or 24 
channels. The energy shift AE of the calculations that opti- 
mizes the fit to the data, and the associated /N , are shown 
in the last two columns. 



R (fm) 


channels 


X 2 /N 


AE (MeV) 


X 2 /N 


3.798 


1 


33.3 


-0.80 


4.65 


3.798 


9 


4.28 


-0.12 


2.71 


3.798 


24 


1.66 


0.0 


1.66 



uncertainty. A possible explanation for the larger radius 
could be the influence of the polarization of high-lying 
states not included in the calculations (see below.) 

It is also encouraging that the extracted values of a r 
shown in Table [TT| are similar to those determined in the 
analysis of the fusion data for 64 Ni+ 64 Ni (a r = 0.403 fm 
3), 16 0+ 16 (a r = 0.41 fm 27]) and 48 Ca+ 96 Zr (a r = 
0.40 fm [H|.) It is noted that in the previous works the 
densities (including the radius) were kept fixed and only 
the value of a r was adjusted in each case to improve the 
fit to the data. 



B. Polarization effects 

The polarization effect discussed in the introduction 
section is illustrated in Table HVl The Table shows that 
for the M3Y+rep-2 potential, one needs to shift the 1 
channel calculation by AE = -0.80 MeV and the 9 chan- 
nel calculation by -0.12 MeV in order to optimize the fit 
to the data. The negative energy shifts are equivalent to 
using a larger radius of the reacting nuclei. For example, 
the required energy shift of -0.12 MeV for the calcula- 
tions with 9 channels can be simulated by increasing the 
radius of 48 Ca by only 0.02 fm. 

The required energy shift AE shown in Table IIVI for 
the calculation with 24 channels is zero simply because 
the radius R was already adjusted in this case to optimize 
the fit to the data. The issue whether the calculations 
have converged with respect to the excitation and po- 
larization of high-lying states is a difficult question to 
answer. It is possible that the polarization of other high- 
lying states, which have not been considered here, could 
play a role and explain part of the 0.05 fm difference 
between the radius of the M3Y+rep-2 solution and the 
estimated experimental matter radius (see Table UlTl ) 



C. S factor representation 

A good way to illustrate the behavior of the fusion 
cross section 07 at low energies is to plot the S factor for 
fusion, 

S = E c . m . 07 exp(r; - 770), (9) 



7 




FIG. 6: (Color online) S factors for the fusion cross sections 
shown in Fig. [5] Also shown is the S factor obtained in 
calculations with the M3Y+rep-l potential. All calculations 
include 24 channels, except the no-coupling limit which has 
only 1 channel. 
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FIG. 7: (Color online) The logarithmic derivative of the en- 
ergy weighted cross sections shown in Fig. [5] Also shown 
is the result of coupled-channels calculations based on the 
M3Y+rep-l potential. The top curve is the constant S factor 
limit, Lcs(E c 



where r\ = Z\Z2e 2 / (Tiv) is the Sommerfeld parameter 
and 770 is that value of rj at a fixed reference energy E . 
The S factors for the fusion cross sections shown in Fig. 
[5] are illustrated in Fig. [5] using the (arbitrary) reference 
energy Eq — 52 MeV. Also shown is the result obtained 
with the M3Y+rep-l potential and 24 coupled channels. 

The coupled-channels calculations for the Woods- 
Saxon potential produce an S factor in Fig. [5] that keeps 
increasing with decreasing energy. The S factors ob- 
tained with the two M3Y+repulsion potentials and 24 
coupled channels are lower. The S factor for best fit to 
the data (the solid curve, obtained with the M3Y+rep-2 
potential) has a maximum at E s = 43.2 MeV. The cross 
section associated with the latter maximum is very small, 
about 0.3 nb. The S factor for the calculation based on 
the M3Y+rep-l potential has a maximum at E s = 35.4 
MeV which is outside the depicted energy range. 

It would be very interesting to know whether the pre- 
dicted S factor maximum near E s = 43.2 MeV can be 
confirmed by experiments but to measure a cross section 
of only 0.3 nb would be a serious challenge. 

D. Logarithmic derivative 

The logarithmic derivative of the energy-weighted 
cross sections, 

L(E c . m .) = -J— ]n(E c . m a f ) , (10) 

is illustrated in Fig. [7J The logarithmic derivatives de- 
rived from the data and from the coupled-channels cal- 
culation based on the M3Y+rcp-2 potential (the solid 



curve) are seen to be in verygood agreement. Similar 
results were obtained in Ref. [5| in coupled-channels cal- 
culations that used a Woods-Saxon potential with the 
large diffuseness a = 0.9 fm. Thus it appears that a large 
diffuseness of the Woods-Saxon has an effect that is sim- 
ilar to that of the shallow M3Y+repulsion potential, at 
least in the low-energy region discussed here. 

The similarity of the M3Y+repulsion and a Woods- 
Saxon potential was recently pointed out in Ref. p9l |. 
It was shown that the M3Y+repulsion potential can be 
reproduced accurately in the barrier region by a Woods- 
Saxon potential with large diffuseness. However, a nu- 
clear potential with a large diffuseness is inconsistent 
with many measurements of elastic and quasielastic scat- 
tering. For example, a recent systematic study of the 
quasielastic scattering of nuclei showed that a realistic 
diffuseness in the range of 0.64 to 0.69 fm is indeed re- 
quired [301 ] . 

It is very interesting to point out that the low- 
energy behavior of the experimental logarithmic deriva- 
tive shown in Fig. [7] is different from the behavior ob- 
served in other systems, in particular in medium-heavy 
systems 0, where the logarithmic derivative usually in- 
creases linearly with decreasing energy and intersects 
with the logarithmic derivative for constant S factor [23[ , 

Lcs{E c . m .) = nv/Ecm.- (H) 

However, there are other systems that exhibit a deviant 
behavior at low energies. For example, the logarithmic 
derivative for the fusion of 36 S+ 48 Ca [3l| also becomes 
rather flat at low energies and it seems unlikely it will 
intersect with the constant S factor limit. In fact, the 
S factor for the fusion of 36 S+ 48 Ca increases slowly and 
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linearly in a logarithmic plot with decreasing energy (see 
Fig. 3 of Ref. [3lj].) 

Another interesting point is that the solid curve in Fig. 
[7] exhibits a maximum near 48 MeV before it rises steeply 
below 44 MeV. It intersects with the constant S factor 
limit at E s — 43.2 MeV, where the associated S factor 
in Fig. [B] develops a maximum. 

VI. CONCLUSIONS 

The fusion data for 48 Ca+ 48 Ca have been analyzed 
using the coupled-channels technique and different ion- 
ion potentials. The analysis based on a standard Woods- 
Saxon potential clearly showed that the data are strongly 
hindered at low energies. By employing and adjusting the 
M3Y+repulsion double-folding potential it was possible 
to achieve an excellent description of the data. 

The best fit to the data was achieved with a nuclear 
radius of 48 Ca that is slightly larger than but still con- 
sistent with the matter radius of 48 Ca. The latter radius 
was determined from the measured rms charge radius and 
the rms neutron radius extracted from an analysis of elas- 
tic proton scattering data. The fact that the extracted 



radius is slightly larger than the matter radius may be 
caused by the polarization of high-lying states that are 
not included in the coupled-channels calculations. 

The entrance channel potential for the best fit to the 
data has a rather shallow pocket, consistent with the find- 
ings of previous analyses of fusion data for medium- heavy 
systems. The M3Y+repulsion potential model is there- 
fore also referred to as the shallow potential model, in 
contrast to models based on the standard Woods-Saxon 
potentials, which have relatively deep pockets in the en- 
trance channel potential. 

The S factor for the fusion of 48 Ca+ 48 Ca does not show 
a maximum within the energy range of the experiment. 
However, it is predicted to develop a maximum at a 3 
MeV lower energy which is nearly the same as the energy 
value obtained from the extrapolation method in Ref. 
[Iol |. The cross section associated with the maximum S 
factor is very small (w 0.3 nb) and is a serious challenge 
to the experimental technology. 
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